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Abstract: Recent measurements of the Belle Collaboration of the exclusive production of two charmonia 
in the e~*"e~ annihilation differ substantially from theoretical predictions. Recently it was suggested that a 
significant part of this discrepancy can be explained by the process e'^e~ — > 2J/il: and value for its cross 
section was published. In this paper it is shown, that these results are incorrect and the values for the cross 
. . , sections of the production of different vector charmonia are presented. 
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^ ■ The charmonium states (such as J/V^, ^/'(25') or r/c) are of great interest both from theoretical and 

^ . experimental points of view because of (i) their clear experimental signature and (ii) great theoretical 

1—5 ! simplifications which arise from their non-relativistic nature. Especially interesting are reactions 

\0 [ of exclusive double charmonium-states production in e"'"e~ annihilations since all non-pertubative 

^ ■ constants can be determined phenomenologically from (cc) — e^e~ decay rates. 
fS| ' Recently the Belle Collaboration has observed e"'"e~ annihilation in two charmonium states at 

^ . center-of-mass energy = 10.6GeV and their measurements differ significantly from the results 

Q\ '. based on non-relativistic quantum chromodynamics (NRQCD) p|, Q . For example the measured 

^ ■ cross section for the process e"'"e~ — > J/iprjc is about one order of magnitude higher then the re- 

^ ■ suits obtained in 0, ^, ^. This presents a challenge to our current understanding of charmonium 

• production. 
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The possible way to explain the part of this difference was proposed in (BBL). The authors 
assumed that some of the events in Belle's J/iprjc signal could actually be J/ipJ/ip events produced 



\^ \ in two virtual photons e"'"e annihilation and obtained the values 

aBBL[e+e- ^ 7/7/- + J/^] = 8.70±2.94fb, (1) 
> : aBBL [e+e~ ^ J/^l) + i){2S)] = 7.22 ± 2.44 fb, 

cTBBLle^e- -^ij{2S)+i>{2S)] = 1.50 ± 0.51 fb. 

H ; 

These results have some drawbacks that follow from the approximations used by the authors: 

(i) The calculations in BBL were held with the help of NRQCD and instead of physical masses 
of vector mesons the values Mj/^ = M^(2S) = 2mc = 2.8 GeV were used. Since in the experiment 
the reaction e^e~ —>■ VV goes near the threshold value of the energy cross section depends strongly 
on final particle masses and this approximation can lead to large errors. 

(ii) In BBL all QCD corrections were neglected. As it was shown in p| the effect of these 
corrections can be significant. For example the leading order value of the probability factor used in 
li^i) j/ip]lo = 0.208 Gev^, but if first order QCD corrections were taken into account it becomes 
[(Oi)^/^]^LO = 0.335 GeVl 

Both of these drawbacks could be corrected if we use not model but physical values of input 
parameters, i.e. vector meson masses and eeV coupling constants. 
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Figure 1: QED diagrams for e^e 27* ^ V1V2 



Four diagrams of the process e~^e~ V1V2 are shown on figure |l|, other two can be obtained 
from diagrams ^.c and l.d by interchanging final state vector mesons. First of all, one can note, that 
diagrams l.a and l.b are enhanced with respect to others by the factor ~ (s/4My)^ ~ 10, which 
arises from virtual-photon propagators. The other source of this enhancement is that diagrams l.a 
and l.b are strongly peaked near the beam direction. The QCD corrections caused by the gluon 
emitted and captured by the same quark or by the quark- ant iquark pair which forms a vector meson 
are taken into account in eeV coupling constant. When gluon line connects quarks of different vector 
mesons we get additional color and photon-propagator suppressions. So from now on only diagrams 
l.a, l.b will be considered. 

The amplitude for this diagrams is 



M 



e 9i92(pi(p2 



+ 



ql 



where 0° is i's meson polarization vector, gi is the coupling constant of the vertex eeVi which can be 
obtained from the electron decay width of the vector meson Vi^2'- 



9^ 



'nee 



Mi is the mass of Vi and F^^ is the widths of its electron decay. 
The differential cross section 
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Table 1: 
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Table 2: (x(e+e- V1V2), fb 
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where x = cos 6 is the cosine of the angle between e~ and Vi momenta, 



l^s - (Ml + M^y^s - (Ml + M^y. 



(2) 



The plot for differential cross section for the case Vi = V2 = J/ip is shown on figure 2. One can 
easily see that, as it was mentioned above, it is peaked near the values x = ±1. 

The result for total cross section is obtained from (|]) by integrating x from -1 to 1: 
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(3) 



In the case of identical final particles (^) should be divided by 2 to avoid double counting. 

Using experimental values of vector meson masses and electron decay widths listed in table 1 
one can obtain total cross-sections presented in table 2. These results differ significantly from those 
presented in BBL, for example for Vi = V2 = J/ip 
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Figure 2: Angular distribution for the process e+e —>■ 2J/ip 

Thus the author would like to note that the explanation proposed in BBL is unsatisfactory and 
the question of the difference between results of the Belle experiment and their theoretical predictions 
remains opened. 

Author would like to thank A.K.Likhoded and S.S.Gershrein for valuable discussions. 
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